New conditional symmetries and exact solutions of 
nonlinear reaction-diffusion-convection equations. Ill 

Roman Cherniha0 and Oleksii Pliukhirll 

Institute of Mathematics, Ukrainian National Academy of Sciences 
Tereshchenkivs'ka Street 3, Kyiv 01601, Ukraine 



A complete description of Q-conditional symmetries of reaction-diffusion-convection 
equation with arbitrary power nonlinearities is finished. It is shown that the results 
obtained in the first and second parts of this work (see arXiv: math-ph/06 12078 and 
arXiv: math-ph/0706.0814) cannot be extended on new power nonlinearities arising in 
the diffusion and convection coefficients. 

1. Introduction. 

In the paper p[] (see [2]- [3] for details) the complete description of Q-conditional symmetries 
of reaction-diffusion-convection equations 

Ut= [[/'"t/.]. + A[/™[/, + C(f/), (1) 

Ut = [U^U^l + XU^+'U, + C{U), (2) 

where A and m are arbitrary constants and C{U) is an arbitrary smooth function, has been 
done. The symmetries obtained for constructing exact solutions of the relevant equations have 
been successfully applied. In the particular case, new exact solutions of nonlinear reaction- 
diffusion-convection (RDC) equations arising in applications have been found. 

The most general RDC equation with power functions arising in the diffusion and convection 
coefficients reads as 

Ut= [U"'u^i + xu^a, + c{u). (3) 

In the cases n = m and n = m + 1, this equation coincides with ([1]) and ([2]), respectively. Here 
we report the main result concerning the structure of Q-conditional symmetries of equation ([3]). 
Note that equation ^ with n = is reducing by local substitution y = x + Xt to the equation 

Q-conditional symmetries of this reaction-diffusion equation were investigated in [4J, so that 
we assume n 7^ 0. 
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2. Main Result. 



This is well-known that a (l+l)-dimensional evolution equation may admit Q-conditional 
symmetries of two different forms 

Q = dt + e(t, X, U)d, + r^{t, X, U)du, (4) 

and 

Q = d^ + r]{t,x,U)du, (5) 

where ^ and r] are unknown smooth functions, which should be found. 

We have proved that equation ([3]) with Xn ^ 0, n m,m + 1 cannot admit any new Q- 
conditional operators of the form (jlj). In other words, a RDC equation with power coefficients 
of diffusion and convection admits Q-conditional symmetries of the form (jl]) only in the cases 
presented in [1]. 

Theorem 1 The RDC equation Ij^ with Xn 0, n ^ m,m + 1 is invariant only with respect 
to operators of the form @), which are equivalent to a linear combination of the Lie symmetry 
operators listed in table 1 of /J/. 

It should be stressed that we didn't consider the problem of constructing Q-conditional 
symmetries of the form ([5]) because one is equivalent (up to the known non-local transforma- 
tion) to solving the given equation ([3]) [6J. Obviously, the nonlinear RDC equation ([3]) is not 
integrable therefore a complete description of Q-conditional symmetries of this form cannot be 
derived. On the other hand, one can try to find particular solutions of the relevant determining 
equation for the function ?], which was derived in [7j, and to construct some operators of the 
form ([5]). 

To prove theorem 1 we use the substitution [2j 

\ln[/, m= -1. ^ ' 

In the case m ^ —1 substitution ([6]) reduces equation ([3]) to the form 

V,^ = VWt~XV''V, + F{V), A^O; A;^0,p,p + 1, (7) 

where p = -j^, k = F{V) = -{m + 1)C {v"^^ , and in the case m = -1 to the 

form 

v.. = e'^Vt - Ae("+i)^14 + F{V), A ^ 0, n ^ 0, -1 (8) 

where F{y) = — C(e^). We use the work [7] to obtain the system of determining equations for 
finding the coefficients of the operator 

Q= dt + at,x,V)d,+r]{t,x,V)dv, (9) 
which is locally equivalent to the operator (jlj) (up to notations). 
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In the case Fo{V) = VP, Fi{V) = -XV\ F2{V) = F{V), system (2.38) [7j takes the form 
Cvv = 0, 



{2^vV - - ^t)VP - ^r]pVP-' - X^.y'^ - \ki^V^-^+ ^0) 
+?>ivF - 2ri^v + U = 0, 

r]Fv + {2^, - Vv)F + {2C,v + VtW^ + prfV^-^- 

-XtJ^V" - T]^^ = 0, 

and in the case Fo{V) = , Fi{V) = -Ae^'^+i)^, F2{V) = F{V), one takes the form 

Cvv = 0, 

= 2ey(-Ae(«+i)^-ee^) + 2e.y, 

(6 + + (e + A(n + 1) - 2^v)v)e''+ (11) 
+Ae.e("+i)^ - S^vF + 2r],v - U = 0, 

r]Fv + (2e. -w)F + ri\^ + 2i,r]e^ + r^^e^- 

We divide the solving systems fllOj) and (ITT]) on three cases: 

(a) ^ = aV + b, r) = ijiV), a = const, b = const 

(b) ^ = a{t, x)V + fit, x), a{t, x) ^ 0, 

(c) i = fit, x), 7] = git, x)V + hit, x). 

One can easily check that these cases take into account all possible solutions of the systems 
(fTOD and (fllD. 



Theorem 2 In the cases (a) and (6) equations ^ and ^ can be invariant only with respect 
to the Lie symmetry operator of the form 

Q = dt + ^d^, C = const. (12) 

Theorem 3 In the case (c) equations ^ and ^ are invariant only with respect to operators 
of the form which are equivalent to the Lie symmetry operators listed in ^ and 

Theorem 1 immediatelly follows from theorems 2 and 3 if one takes into account that the 
RDC equation ([3]) is locally equivalent to the equations ([7]) (if m 7^ — 1) and ([8]) (if m = — 1). 



3. Proof of Theorem 2 



Firstly, let us consider the case (a). In the quite similar way as it was done in [2] (see 
pp. 11-14) one proves that operator (jl]) may take only form f|T2|) . which is the Lie symmetry 
operator of the equations ([7]) and ([8]). 

Case (b). Let us consider the system ffTOj) (the consideration of the system ffTTl) is quite 
similar). The general solution of the first equation arising in flTUl) is 

^ = a{t,x)V + f{t,x), (13) 

the solution of the second equation of ffTOl) is 



n = 2a^ _ 2a/ _ 2Aa T/fc+2 i 



(p+2){p+3) (p+l)(p+2) (fc+l)(fc+2) 

+a^V^ + g{t, x)V + hit, x), 



if p 7^ —2, —3, A; 7^ —1, —2 (the consideration of the cases p = —2, —3, A; = — 1, —2 is much 
more simpler). Substituting (1141) into the third equation of system (llOp we obtain 



F{V) 



1 / _ 2a^(p-2) _ 2pp _ 

3\ (p + 2)(p + 3) {p + l){p + 2) 

_0 \ „ ( ^ L P ~ ^ ] T/p+fc+2 _ 

y{p + 2){p + 3)^ {k + l){k + 2)J 

_2Af + ^ ^ y'^+p+i - 

n(P+l)(P + 2) (A: + l)(fc + 2); 

(p + l)(j9 + 2)(p + 3) p + 2 

+ -f f2(a/, + aj)^^ + at+pfa, + ag{p - 2))^^+^ + 
a V V p + 1 / 

+ (/t + 2//, + pfg + ah{p - 2))VP + pfhVP-' + 
k + 1 

+^''-^ + 2*-^-)-(*TWT2)^"")- <'^) 

Since the function F depends only on the variable V, all coefficients by different powers of this 
variable must be constants. In the general case, we obtain fifteen equations for determining the 
functions a, /, g, h. The number of equations may be shortened for certain values of k and p . 

It turns out that the functions g and h must be constants, if a{t, x) = const, f{t, x) = const 
therefore we arrive at the case (a), which lead only to the Lie operator ([T2|) . Let us prove this 
statement. 

Expression ffTSjl with a = const, f = const, takes the form 



F = UMiV) + -(agip-2)V^^' + ipfg + ahip-2))V^+ 

•J ^ ^ ^ ^ (16) 

+pfhVP'^ + XkgV^ + \khV^~^ + 2g^ 



A 



where 

M(V\ — 2a^(p-2) Ti2p+3 _ 2p/2 T/2p+l 
^l'-' I >^ J - (p+2)(p+3) ^ (p+l){p+2) 

\{p+2){p+3) (fc+l)(fc+2)y 

_9\ f ( A: I p \ T/A:+p+l_ 

l^(p+l){p+2) + (fc+l){fe+2) y 

—An f T/2P+2 _ 2A^fc T/2fc+l 

(p+l)(p+2)(p+3) (fc+l)(fe+2) 

is the polynomial, which depends only on V . Analyzing powers p + 1, — 1, k, k — 1, in 
f|T6|) . we receive the conclusion that there are only five cases, when two and more among them 
are equal, namely: 

k=p-l, k = p + 2, p = 0, p=l, k = l. (17) 

We have to investigate also the general case, when conditions f lT7|) don't take place. In the 
general case, all coefficients by the powers of V must be constant. Let us consider the coefficients 

by the terms V^^^ : -g{p — 2) = const and : — {pfg + CLh{p — 2)) = const. Considering 

3 3a 
separately two subcases, p ^ 2 andp = 2 one easily obtains that g{t, x) = const, h{t, x) = const, 

provided a{t,x) = const, f{t,x) = const. 

Let us consider the case k = p— 1 from ( fTTl) leading to the powers p+l, p, p—1, p— 1, p—2, 
of the variable V in ( |T6l) . We must separately consider the subcases p = 2,p = —l,p = 0,p = 1 
and p 7^ 2,1,0, -1. 

The expression (fT6l) with p = 2, k = p — 1 = 1 gets the form 



F = ^ (m( V) + i (2fgV' + {2fh + \g)V + h + 2g,^ ) . 



'1^ 



Since the right-hand-sight of (fT8|) cannot depend on t, x, we immediately obtain g = const and 
h = const. 

In the case p 2, 1, 0,-1, k = p — 1 expression (fT6|) takes the form 



F = ^(^M{V) + ^(^agip - 2)V^+' + {pfg + ah{p - 2))Vp+ 
+pfhVP-^ + X{p - l)gVP'^ + \{p - l)hVP-^ + 2g^ 



Since coefficients by 1/^+^ and must be constant we again obtain g = const and h = const. 
Subcase p = —1 contradicts to the condition presented above in ([7j), while the subcases p = 
and p = 1 should be considered for arbitrary k (see ( fTTl) ). 

We have checked other cases from ( fTTl) and obtained the same result. Thus, the functions 
g and h are some constants in the expression ( fT5|) provided a = const, f = const. 

On the other hand, the functions a and / are constants if the powers 2p + 3 and 2p+ 1 are 
not equal to any other power in (fT5|) . 

Thus, to prove theorem I2] we must consider only such cases, when the powers 2p + 3 and 
2p + 1 are equal to other power (s) in (fT^ . One can easily check that this happens only in the 
following cases: 

l)p= -4, 2)p= -|, 3)p= -i, 4) p = 0, 5) p = 1, 6) p = 2, 

7) k = p-l, 8) k = p + 2, 9) k = 2p, 10) k = 2p + l, 11) k = 2p + 2, (19) 
12) k = 2p + 3, 13) k = 2p + 4. 



Let us consider the case 1) from (fT^ in details (all other cases can be investigated in a 
quite similar way). Consider fifteen powers of the variable V, arising in fllSI) with k = p — 1: 
2p + 3,2p+ 1, 2p + l,2p,2p + 2,p + 2,p + l,p,p — l,p,p — l,p — 2, 1, 0, 2p — 1. Let us form the 
table 1 with those values of p listed in the second row, when at least one of the powers listed in 
the first row is equal to 2p + 3 (we don't write down the power 2p — l because the corresponding 
coefficient is already constant). 



Table 1. 





2p+l 


2p 


2p + 2 


p + 2 


p + l 


p 


p — 1 


p-2 


1 





2p + 3 








-1 


-2 


-3 


-4 


-5 


-1 


3 
2 



The values p = —1,-2,— 3 contradict to the conditions presented above (see ([7]), (JH])). 
Other values p listed in the second row of table 1 are only subcases of the corresponding cases 
from f|T9|) . namely: p = — | and p = —4 are subcases of cases 7) and 11), respectively, p = —5 
lead to A; = —6, so that this is a subcase of 8). 

Thus, the power 2p + 3 doesn't coincide with any other, so that, a = const (see the first 
term in right-hand-side of f[T^ ). Otherwise we should put p = 2 but this is another case from 

m- 

To establish that / = const we analyze the power 2p + 1. Let us form the table 2 in the 
same way as we built table 1 and taking into account that a = const. We obtain five special 
values of parameter p. 



Table 2. 





2p 


2p + 2 


p+l 


p 


p — 1 


p-2 





2p + l 









-1 


-2 


-3 


1 

2 



Again the values p = 0, —1, —2, —3 contradict to conditions presented above (see ([7]), (|T^). 
The last value P = —\ listed in the second row of table 2 is only a subcases of case 3) from f |T9l) . 
This means that the power 2p + l doesn't coincide with any other power, hence, / = const (see 
the second term in right-hand-side of (ITS!) ). 

Thus, we obtain that a = const and / = const in the case 1) from (fTOl) . so that, taking into 
account the statement proved above, we arrive at the case (a). 

All other cases from 0191) have been analyzed a similar way and the same result established, 
i.e. a = const and / = const. 

The proof is now completed. 

4. Sketch of the Proof of Theorem 3 

Substituting ^ and r] from the case (c) into the third equation of ( fTOl) (the investigation of 
system ( |TT|) is analogous), we obtain 

(2//. + ft+ pfgW + pfhV''-^ + A(/, + kg)V'' + XkhV^-^ + 2g, - = 0. (20) 

To analyze fl20|) one needs to consider only the following special cases: p = 0, p = 1, k = 
p — 1, k = 1 and the general case, when p and k don't satisfy these restrictions. 



Consider the case p = in detail (the next two cases and the general case are investigated 
in the same way). Equation (l2Up with p = takes the form 

A(/. + kg)V'' + XkhV''-' + ft + 2fU - f.. + = 0. (21) 

Since k p = and k ^ p + 1 = 1, one can split (!2T|) by different powers of V and obtain the 
system 

Xkh = 0, 

A(/. + kg) = 0, (22) 
ft + 2//, - + 2g., = 0. 

Taking into account the restrictions presented above, we obtain 

h = Q,h = -kg. (23) 
Substituting expressions (!23|) into the fourth equation of (fTOj) . we arrive at the linear ODE 

gVFv - {2k + l)gF = Xg.,V'+' + (fc + 2kg^ - gt)V (24) 
with the general solution 

r2k+l ^9x^^k+l t f 9t — 9xx 



F = Ai\/""+^-^y"+^+ r^ " -g)V. (25) 

Note the special subcase (7 = immediately leads to /i = and / = const, therefore the Lie 
symmetry f|T2|) is obtained. 

Since right-hand-side of (l25ll cannot depend on the independent variables, we obtain 

9t-9xx _ „ _ \ y^"^) 

where lambda-s are some constants. The general solution of fl26|) is g = a{t) and then, using 
fl23|) . we have 

/ = -A;a(t)x + /3(t), (27) 

where a()f:), /5(t) are to-be-determined functions. Substituting flTTI) into the third equation of 
we obtain 



{-kat + 2A;2a2)x + /3t - 2A;a/3 = 0. (28) 
Obviously, equation (!28|) is equivalent to two ODE equations with the general solution 
1 A2 

a = ; — , 3 = — — , A = const, i = 1,2. (29) 

2kt + Ai' ^ 2kt + Ai' ' ' ' ^ ^ 

Finally, taking into account (1231) . (l25ll . ( 1271) and ( l29l) . we have found that equation 
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is invariant under the operator 

Q = ^^ + 01,] A ((^"^ + ^2)5x - m 

2kt + Ai\ 

Multiplying this operator by 2kt + Ai one obtains the operator 

X = {2kt + Ai)dt + {kx + A2)d^ - Vdv, (30) 

which is nothing else but a linear combination of Lie symmetry operators of this equation (see 
the case 7 of table 1 [5]). 

Let us consider the case k = 1, which is special. Equation fl20l) with A; = 1 takes the form 

(2//.. + ft+ pf9)V' + pfhV^'' + A(/,. + g)V + Xh + 2g, - = 0. (31) 

We must consider only two subcases p 2 and p = 2 (we remind the reader that the subcases 
p = and p = 1 were considered above). It is easily shown that the first of them leads only to 
Lie symmetries, since we can split the equation fl3Tl) with respect to the four different powers 
of V, i.e., four equations are obtained. In the second subcase ( !3T|) takes the form 



(2//. + ft + 2fg)V^ + {2fh + A(/,, + g))V + Xh + 2g., - = (32) 
and we obtain only three equations 
2//. + ft + 2fg = 0, 

2fh + X{f, + g) = 0, (33) 
Xh + 2g^ - f^^ = 0. 

So, (l32l) is the nonlinear system of three PDEs for three unknown functions. 

We must also take into account the fourth equation of (fTOj) with p = 2, namely: 



{gV + h)Fv + (2/, - g)F = -{gV + h){?>gV'' + 2hV)- 

-(2/. - g){gV + h)V^ - {g^V + ht)V^ + X{g^^ + h^)V + g^^V + h,^. 



(34) 



Analyzing the first term of fl34l) we must consider three subcases, 1) (7 = /i = 0, 2) = 0, h ^ 
0, ?>) g ^ 0, for solving this equation. It turns out, subcases 1) and 2) lead the Lie symmetry 
f lT^ (we omit the relevant calculations because they are rather simple). 

Let us consider the most complicated case 3) 51 7^ 0. The general solution of the equation 
([31 is 



Agl^^ + AsF^ + Ai"\/ + Ao + A (^y + \ A = const. (35) 



Here Aj, i = 0,1,2,3 are some functions on /, g, h and their derivatives but we omit those 
expressions for lambda-s because they are rather cumbersome. Since the function F must 
depend only on V, then we obtain 

1 - ^ = A* A* = const (36) 
9 
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if A 7^ 0. Formula (1361) may be written in the form 

g = A,U,, Ai = ^^ (37) 

(the case A\ = 1 leads to the contradiction g = Q). Substituting ( 1371) into the third equation 
of (133|) we obtain 

h = ^^^/... (38) 

Substituting (1571) and (1551) into the second equation of (1551) . we arrive at 

2(2Ai - 1)//,, = \\A^ + 1)/,. (39) 
Equation (139|) is reduced by the substitution y(/) = /^^ to the form 

where 7(t) is an arbitrary smooth function (the case Ai = ^ leads again to the contradiction 
g = 0). Substituting (HUj) into the first equation of (155]) and taking into account (1571) . we obtain 

Differentiating (HOj) by t and (HTl) by x and equaling the expressions obtained, we have 
m + l)/7.^AV.H-iP^, A^^^ (42) 



2^1 - 1 V27 '2{2Ai - 1) 
The expression ( H21) is satisfying only with Ai = —1, but this leads to = (see (HU]) ) and 
therefore 5^ = (see (137|) ). Thus, the contradiction is again obtained. 

Let us consider (l35l) with ^4 = 0. Substituting (135|) with A = 0, into equation (15^ and 
splitting expression obtained by the different powers of V, we obtain the system 

gt + 2ig + X;)ig + f,,) = 0, 

2Uh + A2) + giAh + A2) + 3Xsh + - Xg^ = 0, 

2/i(/i + A2) + 2Ai/,.-A/i,-fc = 0, ^^^^ 
Ai/i + 2Ao/x - Xog - h^^ = 0. 

Finally, any solution of the nonlinear system (l33l) and (H3|) . which consist of seven equations 
on three functions /, g and h generate the operator of Q-conditional symmetry 

Q = dt + fd, + {gV + h)dv (44) 
for the equation 

Ku. = VWt - XVV^ + XsV^ + X2V^ + XiV + Ao. (45) 

It turns out, that the overdetermined system of PDEs (155]) and (H5]) is compactable. How- 
ever, all its solutions produce the operators of the form (IH|) . which are nothing else but Lie 
symmetry operators of (H5l) obtained in [7]. We have established this using computer algebra 
package Mathematica 5.0. The relevant calculations are omitted because their awkwardness. 

The sketch of the proof of theorem 2 is now completed. 
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